Entanglement for rank-2 mixed states 
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In a recent paper, Rungta et. al. [Phys. Rev. A, 64, 042315, 2001] introduced a measure of mixed- 
state entanglement called the I -concurrence for arbitrary pairs of qudits. We find an exact formula 
for an entanglement measure closely related to the /-concurrence, the I-tangle, for all mixed states 
of two qudits having no more than two nonzero eigenvalues. We use this formula to provide a tight 
upper bound for the entanglement of formation for rank-2 mixed states of a qubit and a qudit. 
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Quantum entanglement, as typified by the singlet state 
l^ - ) = "^7|(|01) — 1 10)), is a uniquely quantum property 
of any nonlocal superposition state. Unlike other proper- 
ties of quantum states, entanglement is singled out by the 
apparently central role it plays in quantum information 
processing. Indeed, it is becoming clear that entangle- 
ment is a resource which can be used to perform quan- 
tum communication and computational tasks jlj] . To this 
end, it is of paramount importance that a general theory 
of entanglement be developed in order to understand how 
it can and cannot be manipulated. The first step in this 
program is to develop physically motivated entanglement 
measures to quantify how much entanglement is present 
in a given state. Many entanglement measures have been 
proposed, however, for mixed states, they typically in- 
volve difficult minimisations. 

The development of an efficient formula for the entan- 
glement present in arbitrary mixed states is a crucial goal 
of quantum information science. Such a formula would 
lead to great advances in the theory of quantum com- 
munication, as well as possibly revealing the correct role 
entanglement plays in quantum computational speed-up. 
To date, while there has been a great deal of work on 
this subject, there is only one easily-computable entan- 
glement measure and it is specific to pairs of qubits ||, |) . 
The purpose of this Letter is to add to this situation by 
providing a formula for one particular bipartite entangle- 
ment measure, the /-tangle, for mixed states which have 
at most two nonzero eigenvalues. 

The structure of this Letter is as follows. We begin by 
introducing the concurrence and the tangle for a pair of 
qubits. We then discuss a generalisation of the concur- 
rence and the tangle for pairs of qudits. The main result 
of this Letter, a formula for the /-tangle of rank-2 mixed 
states, is then established. We conclude with a corollary 
of the main result, an upper bound for the entanglement 
of formation of a rank-2 mixed state of a qubit and a 
qudit. 

Before we discuss the /-tangle, we introduce the con- 
currence, a mixed-state entanglement measure for states 
of a pair of qubits AB g, |, £§. The definition of the 
concurrence makes use of a specific transformation on 



density operators, the spin-flip operation, which is de- 
fined as follows. Consider an arbitrary mixed state p of 
AB. We define the spin- flip of p to be 



p = tr(p f )/«) / - p\ 



<Pb 



(1) 



where pa — ^b(p) and ps = ^a(p) denote the reduced 
density operators for subsystems A and B respectively. 
(We have included the trace and hermitian adjoint terms 
so that the spin-flip operation is defined for arbitrary op- 
erators acting on AB.) The formula for the spin-flip is 
applicable to arbitrary bipartite systems, in which case 
it is called the universal state inverter ||. 

The spin-flip operation Eq. ([!]) on a pair of qubits is an 
example of an antilinear operation. To be more precise, 
consider a pure state p — (ip\. The spin-flip operation, 
when japplied to this state, is equivalent to the expression 
p = \ip)(ip\, where 



\i/>)=<r v ®<r v {\rl>))*, 



(2) 



and where o~ y is expressed in the computational basis 
as ( 9 ) , and the complex conjugation is taken in the 
computational basis. The operation in Eq. (Q) is clearly 
an antilinear operator The definition of the spin- 

flip as an antilinear operator extends, via linearity, to 

all mixed states, p = OpO, where we have added the 
arrows above the antilinear operator 6 representing the 
spin-flip to indicate the direction in which it acts. It is 
worth noting that the description of the spin-flip Eq. (Q) 
in terms of an antilinear operator 9 is specific to two 
qubits. 

For pure states p = \ip)(ip\, the concurrence C of 
is defined to be C = = (ip\p\ip) . When the 

state p of the two qubits is mixed, the concurrence C is 
defined to be a minimum over all pure-state decomposi- 
tions {pi,\tpi}} of p: 



C(P) = r min y2pi\{i/>i\il>i 



(3) 



It is convenient to introduce another entanglement mea- 
sure closely related to the concurrence, the tangle r [n2| , 
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which is also defined as a minimisation over pure-state 
decompositions: 



r(p) 



mm 

{pi,\M>i)} 



Pi\&i\Tpi)f 



(4) 



The squared-concurrence satisfies the inequality C 2 < 
t, which follows from the convexity of | (V'il^/'i) | 2 = 
C 2 (|-0i))- It turns out that the reverse inequality also 
holds, so that the tangle is equal to the square of the 
concurrence, r( ( o) = C 2 {p) (The reverse inequal- 

ity may be established by noting that there exists a de- 
composition {pi, achieving the minimum in Eq. (j^) 
which has the property that C(\ipi)) = CQipj)) §]. The 
inequality follows from substituting this decomposition 
into the expressions for r and C 2 .) A simple formula for 
the concurrence of two qubits is known [pi, 



C(p) = max[0, Ai - A2 - A3 - A 4 ], 



(5) 



where the A^ are the square-roots of the singular values, 
in decreasing order, of the matrix pp. 

For a pair of qudits AB we use a variant of the I- 
concurrence of Rungta et. al. Q to measure the entan- 
glement for mixed states of A and B. The /-concurrence 
is defined via Eq. (0) f|, 



C(p) 



min > 

{p<>IV-i>}4^ 



(6) 



where pt = \ipi)(ipi\- The entanglement measure we use 
is a generalisation of the tangle, the I-tangle, defined by 



min y^Pi(ipi\pi\ipi 

{Pi,\lP,)} ; 



(J) 



The /-concurrence and the /-tangle are good mixed-state 
entanglement measures because they satisfy the standard 
properties usually regarded as essential for a good entan- 
glement measure (see, for example, HQ). The inequal- 
ity C 2 < r may be established, by convexity, as for two 
qubits. Because the equal-entanglement decomposition 
only exists for pairs of qubits, the /-tangle is not, in gen- 
eral, equal to the square of the /-concurrence. Based 
on the results of this paper, and the calculations of the 
/-concurrence for isotropic states [[l4| we feel that the 
/-tangle, as defined by a minimisation, is the proper gen- 
eralisation of the tangle Eq. (||) . 

The universal state inverter Eq. (]l|) may be expressed 
in terms of another formula which will be most useful in 
the following. Before we write down this formula, how- 
ever, we need to introduce some definitions. Let \i)a and 
\j)b denote the computational basis states for subsys- 
tems A and B, with dimensions d A and ds, respectively. 
For an arbitrary pair {\i) A ,\i') A }, {\j) B , \?)b} of the 
computational basis states of AB we set up the projec- 



tors P 



(«') 



liMil+li'Mi'l,^ 



Uf) 



\j)B(j\+\f) B (n 



and Q a = P A " ^ ®Pg \ where a = Consider 



the object p a = Q a pQ a - The operator p a is a positive 
operator supported on a 2 x 2 subspace of the Hilbert 
space of AB spanned by {\ij), \if), in this 

way we can think of p a as a subnormalised state of two 
qubits. The two-qubit spin flip, when applied to p a , gives 

Pa = ^apOa = <T V <8> (J V {Q a pQ a)* & V <8> o\ (8) 

where 9 a = 9Q a is the antilinear operator representing 
the spin-flip operation on the 2x2 subspace, and a y is 
naturally defined on the two-dimensional subspaces of A 
and B respectively. Using these definitions we can write 
an alternative formula for the universal state inverter, 



P 



E 



'apt 



(9) 



where the sum over a runs over all of the (d A {d A — 
l)/2)(rfs(ds — l)/2) possible choices of pairs of compu- 
tational basis states. (The reader may verify that Eq. (^) 
follows from the expression of the universal state inverter 
as a tensor product of two superoperators of the form 
V o T. See || for further details.) 

It is convenient, at this point, to introduce two quanti- 
ties which will simplify the statement of our main result. 
Let p be a density operator for a pair of qudits having 
no more than two non-zero eigenvalues. We may write p 
in terms of its eigenvectors, 

P = p\v 1 )(v 1 \ + (l-p)\v 2 )(v 2 \. (10) 

Using these eigenvectors we construct the tensor 

Tijki = tr(7y7 fci ), (11) 

where 7^ = |i>i)(fj|. We also construct the real symmet- 
ric 3 x 3 matrix My whose independent entries are given 
by 
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Mil = 


-Z1221 + 


-Z1122 + 


JT2112, 


M 12 = 


I il221 ~ 


7/2112, 












M13 = 


-T n2 l ~ 


^■/2122 + 


7/1112 — 7/1222 






4 4 










M 22 = 


-4T1221 


+ 2^1122 


— ^-/2112, 








i i 


M 23 = 


4 J n2i- 


4 Tl112 + 


7/2122 — 7/1222 










M 33 = 


4 iim ~ 


2^1122 + 


-J- 2222- 



(12) 



(The entries of M will be shown to be real in the follow- 
ing.) 

We now have all the necessary ingredients required for 
the statement of our main result. 

Theorem 1. Let p be any density operator for a pair AB 
of qudits, of dimensions d a and ds, respectively, having 
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no more than two nonzero eigenvalues. The I-tangle t 
between A and B is given by the expression 

t( P ) = tr(pp) + 2A min (l - tr(p 2 )), (13) 

where A m in is the smallest eigenvalue of the matrix M 
defined by Eq. (]Z|). 

The formula Eq. (|l3|) for the /-tangle is easy to com- 
pute for all rank-2 mixed states of a pair of qudits. We 
have verified this formula by performing the minimisa- 
tion in Eq. (^) numerically for a large number of ran- 
dom states, and compared the results with the formula 
Eq. (|l~3|). We found agreement to within 10~ 16 . 

Proof. The method we use to prove this theorem is sim- 
ilar to that employed by Hill and Wootters 0. 

Consider an arbitrary pure state which can be writ- 
ten as a linear combination of the two eigenvectors of p, 
\ijj) = ci\vi) + c 2 \v 2 ). The /-tangle of \if>) is given by the 
expression 

t(iI>) = (4>\a\ij) = J2^\6 a W(ip\*0 a W, (14) 

a 

where a = \ip)(ip\, and we have used Eq. (0) to rewrite 
the spin flip in terms of the antilinear operators 9 a . Each 
of the terms in the sum over a may be written as a trace, 

tW>) = 5>(^C q ^C q *), (15) 

a 

where u>ij — etc* is the density matrix of l^) expressed 

in the \v 2 )} basis, and Q = (vAQ a \vj). 

The function on the RHS of Eq. ( |15| ) can be extended 
via linearity to a function / of all 2 x 2 density matrices 
lo expressed in terms of the {|fi), \v 2 )} basis, i.e., f(to) = 
J2 a * r {u* (, a Lo(, a *). The function / has the property that 
it is equal to the /-tangle for all pure states tp, f(ip) — 

Any 2x2 density operator u> may be expressed in 
terms of the Pauli matrices via the operator expansion, 
u> = i(Z + r • <x), where = tr(u>a l ). Substituting this 
expansion into the expression for / gives the quadratic 
form 

/(w) = J tr(T) + J2 rjLj + ]T r,r k M jk , (16) 

3 3,k 

where T = ^ a C Q *C Q , 

Li =tiVT), (17) 

and 

Each of the terms in the sum over a in Eq. ( |l8| ) is a real 
symmetric matrix, so that M is a real symmetric matrix. 



It may be straightforwardly verified that the entries of 
M are given by Eq. (p"2|). 

For the rank-2 density operator p, the state-space of 
the system AB can be considered to be the space of all 
convex combinations of superpositions of |i>i) and \v 2 ). 
If a particular state \ip) of AB is pure, its corresponding 
2x2 density operator in the \v 2 )} basis, cu = \{I + 

r • <t), satisfies the condition |r| 2 = 1. In this way, we can 
think of the entire state-space as the Bloch sphere where 
the poles are the eigenvectors \v\) and \v 2 ). A particular 
decomposition of p may be viewed as the weighted sum 
of points on the surface of the Bloch sphere, where p lies 
at the centre of mass of the weighted sum. The function 
/ is defined on the entire state-space |r| 2 < 1. 

When the bipartite system AB is a pair of qubits, there 
is only one term in the sum Eq. (||), and / reduces to 
the quadratic form that Hill and Wootters Q study. In 
this case, the eigenvalues of M are given by ±i|det£| 
and ^tr(C*C), which means that / is convex along two 
directions and concave along a third. In general, the 
matrix M will have three positive eigenvalues, so that / 
is typically convex. 

For the purposes of this proof it is essential that a 
quadratic form g be constructed which agrees with / on 
pure states "0 which has the additional property that it is 
convex along two directions and linear along a third. A 
function g which has these properties may be constructed 
from / as follows, 

5 H^/H-A min (|r| 2 -1), (19) 

where A m i n is the smallest eigenvalue of the matrix M . 
This function is a quadratic form, 

g(u>) = K + Y, r 3 L, + rjr k N Jk , (20) 

3 j,k 

where N = M - A min /, and K = |tr(T) + A min . The 
matrix N which defines the quadratic form g has two 
positive eigenvalues and one zero eigenvalue so that g is 
convex along two directions and linear along the third. 
The quadratic form g has the additional property that it 
is equal to / for pure states ip, (|r| 2 = 1). 

At this point we recall a theorem due to Uhlmann [||, 
which concerns functions of density matrices expressed 
as minimisations over all pure-state decompositions. 

Theorem 2. Let G be a positive real- valued function de- 
fined on pure states. The function S^ , defined for all 
mixed states p, given by 

Sf(p) = min 5>G(Vi), (21) 

where the minimisation runs over all pure-state decom- 
positions of p, {pi, IV')}; is the largest convex function 
which agrees with G on pure states p = ip. 
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The /-tangle is expressed as a minimisation over all 
pure-state decompositions of a density operator, so if we 
could find the largest convex function which agrees with 
/ = t on pure states ?/> it is guaranteed to be equal to the 
/-tangle for all mixed states. We claim that g is precisely 
this function. To prove this statement we proceed via 
contradiction. Assume that there is a convex function g' 
which agrees with / on pure states but which is larger 
than g. We can write g' in terms of g and a 'correction' 
t 

gf(u)=g(u)+t(u>). (22) 

In order that g' satisfy the correct boundary conditions 
t(ui) must be zero on the set of all pure states (the surface 
of the Bloch sphere). Consider the coordinate system Xi 
defined by the three eigenvectors of TV where, without loss 
of generality, we choose the x$ coordinate to correspond 
with the direction along which g grows linearly. The 
condition that g' be convex is equivalent to requiring that 
the Hessian matrix 

H = (23) 

OXiOXj 

is positive semidefinite ||. In order that H be posi- 
tive semidefinite, it is necessary (although not sufficient) 
that the entry H 33 satisfies the inequality H 33 > JlO| . 
Because g grows linearly along z 3 , the condition that 
H 33 > becomes > 0, which means that t is convex 

as a function of x 3 . The boundary conditions, t(ip) = 0, 
therefore imply that t is everywhere negative. Hence 
g' < g, and we have a contradiction. This implies that 
g is the largest convex function which takes the values 
t(V0 on the set of all pure states. Therefore, theorem || 
shows that g is equal to the /-tangle. 

The expression for the /-tangle r may be simplified by 
noting that, for rank-2 p, 1 - |r| 2 = 2(1 - tr(p 2 )). Note, 
also, that f(p) — tr(pp). Hence we can write r(p) = 
tr(pp) + 2A min (l - tr(p 2 )). □ 

The decomposition which achieves the minimum for 
the /-tangle Eq. (|l^) consists of two terms. In contradis- 
tinction to the case of two qubits, the minimising de- 
composition will, in general, consist of terms with dif- 
fering values of r. This is because the surfaces of con- 
stant g will typically be curved, so that the trick of Hill 
and Wootters cannot be applied (see Q for the construc- 
tion of the minimising decomposition when the surfaces 
of constant g are elliptic cylinders). The construction 
of the minimising decomposition follows from observing 
that the function g has the property that it grows linearly 
in one direction. Consider the line parallel to the eigen- 
vector of N, whose associated eigenvalue is zero, which 
passes through the density operator p. The density op- 
erator p may be written as a convex sum of the two pure 
states \cj>i) and \<p2) which lie at either end of the line, 



p = qi\(px) {<f>i\ + ^2102) (02 1- Because the /-tangle is con- 
vex, we obtain the inequality 

t(p) < gir(<£i) + g 2 T((fc). (24) 

However, t = g varies linearly in this direction, so that 
the inequality in Eq. (pi| ) is actually an equality. 

When one of the subsystems of the bipartite system 
AB is a qubit it is possible to obtain a relation between 
the /-tangle r and the entanglement of formation J?. 
For pure states ip of AB the entanglement of formation 
is given in terms of the /-tangle via 

= *(t(i/0), (25) 

where S"{x) = H (A + ^Vl ~ x ) 1 an d H is the binary 
entropy function H(x) — —a; log a; — (1 — x) log(l — x), 
where the logarithm is taken to base 2. The function 
§ is concave and monotone increasing. If we consider 
the minimising decomposition {qi, \(f>i)} we constructed 
in the previous paragraph, we obtain the chain of in- 
equalities 

&{p) < «i*(t(&)) + g 2 ^(r(0 2 )) < *(t(j>)), (26) 

where the first inequality follows from the definition of 
the entanglement of formation, and the second from the 
fact d>(g) is concave along the line passing through the 
pure states \<f>%). This statement is the content of the 
following corollary: 

Corollary 1. For rank-2 mixed states p of a qubit A and 
a qudit B, the entanglement of formation ^ of p satisfies 
the inequality 



^(p)<H(- + - V I Z ^(p)J- (27) 

Numerical experiments indicate that the expressions 
on the LHS and RHS of Eq. @ usually differ only by 
about 10 -4 , so that the inequality is typically very close 
to an equality; it is not, however, an equality. 

We have found a formula for the /-tangle for all states 
of a pair of qudits having no more than two nonzero eigen- 
values. We have also found an upper bound for the en- 
tanglement of formation between rank-2 mixed states of 
a qubit and a qudit. The method we employed to con- 
struct these formulae relied crucially on theorem |^. It is 
tempting to conjecture that the construction of entangle- 
ment measures might be simplified by searching for large 
convex functions which satisfy certain boundary condi- 
tions — perhaps as solutions to certain partial differen- 
tial equations. There are some difficulties with this sug- 
gestion, however. We must include in this search space 
functions which are not analytic at one or more points 
(the concurrence is not analytic at all points). Nonethe- 
less, perhaps there exists a systematic way to construct 
good analytic approximations to entanglement measures 
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expressed as minimisations over pure state decomposi- 
tions. 
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